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Abstract. 

We calculate the Hartree-Fock energy of a density-wave in a spin polarized two- 
dimensional electron gas using a short-range repulsive interaction. We find that the 
stable ground state for a short-range potential is always either the paramagnetic state 
or the uniform ferromagnetic state. The energy of a density-wave state is, however, 
reduced by a factor proportional to (1 — C^), where C is the polarization of the 
electron gas. Since this situation occurs in the most unfavorable conditions (short- 
range repulsive interaction) it is therefore conceivable that by including higher order 
many-body corrections to the interaction a density-wave ground state is indeed found 
to be stable. 
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The possible existence of a non-uniform ground state for a two-dimensional (2D) 
electron system in the presence of a magnetic field has been proposed in explaining 
the results of a number of recent experiments. Several examples of such puzzling data 
are the unidirectional charge- density- wave states that appear in partially filled Landau 
levels m E], the incompressible, inhomogeneous insulating phase in p-GaAs/AlGaAs 
heterostructures displaying metal-insulator transition OHj, and the unusual magneto- 
optical properties of the ferromagnetic phase of p-type Cdi^xMrixTe quantum wells 

Numerical calculations [H^ have definitely proved that at low densities an unpolarized 
electron gas forms a Wigner crystal. However, a study of point defects in the two- 
dimensional Wigner crystal suggests that the quantum melting could be continuous 
rather than first order [Jj, leaving open the possibility of inhomogeneous intermediate 
phases. 

The idea of a non-uniform ground state in a paramagnetic 2D electron system 
has surfaced before. Prior to the discovery of the quantum Hall effect, it was argued 
that in 2D GaAs-type structures the ground state was a charge-density- wave |S]. More 
recent experimental^ E] and theoretical [HI CS] results point indeed to the existence 
of charge- density- wave states in partially filled higher Landau levels on account of the 
quasi-one-dimensional electron motion. 

In this work we are interested in the possible formation of density waves in a 2D 
electron system that exhibits a very large Zeeman splitting. This situation occurs in 
II- VI dilute-magnetic semiconductor structures where the effective Lande factor is up 
to thousands of times the band value. This is exactly opposite to the case in GaAs, 
where the cyclotron energy is dominant. For example, in GaAs the ratio of the cyclotron 
frequency to the Zeeman splitting is around 14, while the same ratio is close to 1/20 in 
Zui^^Cd^Se. 

Quantum Monte Carlo results show that the exchange-correlation hole is larger in 
the polarized electron gas[TTl [T2] suggesting that a highly polarized system is prone 
to density instabilities. Whether or not the instability develops before the Wigner 
crystallization and its dependence on the degree of polarization has not been studied 
yet. Here, we calculate the Hartree-Fock (HF) energy of a spin density wave state in 
a spin polarized system for a delta-function repulsive interaction, the most unfavorable 
situation for the development of a density-wave instability in an isotropic system. 

It is well known that in three dimensions, within the HF approximation and for 
an unscreened Coulomb interaction, the paramagnetic state is unstable with respect 
to formation of a spin density wave with wave vector near 2kF ^H|- However, for a 
short-range potential the stable HF solution is always either the usual paramagnetic 
state or the uniform ferromagnetic state ^3]. We find that this result holds also in the 
2D polarized system. In this case, however, the difference in energy between a density 
wave state with momentum q, E{q)/N, and the ferromagnetic state, i?(q = 0))/A^, is 
reduced by a factor proportional to (1 — C^), where ( = (n^ — ni)/n is the polarization 
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of the electron gas, 
E(q) 



q 



(1) 



If the background charge is allowed to relax, as in the deformable jellium model, our 
result applies to both spin- and charge-density waves 1 1 5j] . 

Since the effective electronic interaction in real systems is somewhere in between 
the short-range delta-function and the unscreened Coulomb potential our result suggest 
the possibility of an inhomogeneous density wave ground state in highly polarized two- 
dimensional electron systems. If this inhomogeneous state exists it will be most probably 
a charge density wave, since electronic correlations in real systems favor charge- density 
over spin-density wave instabilities [^. Below we present the details of our analysis. 

A density wave develops when the correlation function between an electron with 
momentum k and spin a and another electron with momentum k-|-q and spin a', 
(V'ko-V'k+qo-') becomes finite. When a = cr', a charge- density wave is formed, while 
a a' corresponds to a spin-density wave. The interacting electron system Hamiltonian 
is diagonalized by a canonical transformation that introduces a new set of operators: 



lower 



^upper 



COS (6'k/2)^/'k-a<7 + sin (6'k/2)V^k+aa' 
- sin (^k/2)V^k-aa + cos (^k/2)V^k+§<7' 



(2) 
(3) 



where 6'k is the coupling parameter and ip^^'^'^^ {ip^^^'^'') refers to the new lower (upper) 
band excitations. 

The Hartree-Fock ground state energy is a function of q and the parameters ^k 
[Hj . For a system with A^ electrons, volume V and polarization ( the total ground state 
energy is: 



k^ + :^ - (k-q) cos^k 



Eid) =—y 

k 

-^E^(k-k')cos2[(^^k-^k')/2] 



R 



^cos^k- A^C 



(4) 



where the first term is the kinetic energy, the second is the effective Zeeman energy:}:, 
needed to fulfill the constraint of constant polarization, and the last one the Coulomb 
and exchange energies. 

For a delta-function interaction with strength f (k — k') = vq the total ground state 
energy in atomic units § becomes: 




(k • q) cos ^^k 



sm fk 



H 



^cos^k - A^C 



(5) 



\ The spin magnetization is written as M — = — cos^k- 

k 

§ k and q are in units of inverse Bohr radius {(Iq^) and the energy in rydbergs. 
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Figure 1. Normalized Hartree-Fock energy of the spin density wave state, {E{q) — 
E^^^^) / [N Ep) versus its normalized momentum, q/kp, for ( = 0.5 and different 
values of the effective potential strength, We//, as indicated in the legend. 



E{q) reaches a minimum with respect to when 



tan 9h 



2k-q + H'{vo,0 



and 



(6) 



where n = N/V, ( 



AT ^ 



cos 9k and H' = H + Vo(n. 



The quasi-particle energy in the ground state when only the lower band is occupied 
is given by: 

(7) 



lower 



4 2 



The solution for q' = Aq, v'q = vq/ \,V' = V/ is homologous to that for q, vq, V, 
and results in k' = A^/t and H' = X^H. Then, from a numerical point of view, the 
simpler way to minimize the total energy is to take h{vq, () = 1 and a fixed value for 
the parameter H', select the maximum occupied energy in the lower band, and then 
compute n, (, Vq and E{q) for this choice of parameters using Eqs. (0) and 

Fig. ^displays our results for the normalized energy per particle of the spin density 
wave state, E{c\)/{NEf), as function of normalized momentum, q/kp, for ( = 0.5 and 

several values of the effective potential strength, Veff = — = — . The energy is 

Ep 2n 

measured with respect to the energy of the paramagnetic state with the same density 

j^PARA ]^ 



and polarization. 



^,1 + C^) + - C^)- For values of v^jf > 2 there is 

N Ep 2 4 

a range of momenta where the density wave states are lowest in energy than the non- 
magnetic state. However, always the absolute minimum is at q = 0, the ferromagnetic 
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Figure 2. Normalized Hartree-Fock energy of the spin density wave state, {E{q) — 
E^^^-^) / [N Ep) versus its normalized momentum, q/kp, for an effective interaction 
strength of Veff = 2.15 and C = (stars), ( = 0.5 (circles) and C = 1 (squares). 



state. The range of momenta where the density wave states are stable increases with 
the strength of the interactions, while the energies get reduced. 

Fig. 121 shows the normalized energy for Veff = 2.15 and three values of the 
polarization, ^ = 0, 0.5, 0.9. The range of momenta where the density wave state is 
lower in energy is the same for all values of polarization. However, the energies increase 
with polarization. 

By scaling the energies and momenta we find that all our numerical results overlap 
on a single curve. Fig. (jH]) displays our results for the energy per particle relative to 
the energy of the paramagnetic state Ae(q) = (-E(q) — E^^^^)/{NEp) divided over 

^sci = {veff — 2)(1 — C^)/4 as a function of the scaled momentum q = . We 

kp^Veff - 2 

have studied a full range of values of vq and (. All the data are in excellent agreement 
with the surprising simple scaling relation: 

^ = g--l. (8) 

Therefore, the HF energy of a two-dimensional density wave (Eq. (jSJ) is a quadratic 
function of the wave- vector for any value of the polarization and strength interaction. 

In conclusion, we have found that independently of the polarization, for Vejf > 2 
and wave- vector q < kp^yvejf — 2 the density wave ground state is lower in energy 
than the paramagnetic state. In all cases, the absolute minima is the ferromagnetic 
state with q = 0, but the depth of the minima scales with the polarization as (1 — C^). 
As a consequence, in a highly polarized electron gas the density wave states are very 
close in energy to the ferromagnetic instability. Because a short range potential is the 
most unfavorable for the development of density instabilities, that fact suggests that in 
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Figure 3. Scaled Hartree-Fock energy per particle of the spin density wave state, 
Ae{q)/€sci, as a function of the scaled momentum, q = q/{kp^Veff — 2). Circles 
correspond to different values of Veff (between 2 and 8) and polarization C — 0,0.5 
and 0.9. The solid line corresponds to the scaling relation Eq. (jS)). 

real systems the ground state might be a density wave state in some range of densities 
and spin polarizations I6|. Since in two dimensions transitions to long-range order 
phases take place at zero temperature, the extension of our arguments to the geometry 
of real systems, quasi-two-dimensional heterostructures, thus allowing direct comparison 
with experimental results [H], is currently being investigated. 
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